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Abstract 

We discuss the concept of characteristic squeezing modes applied to a travelling-wave optical 
parametric amplifier pumped by an ultrashort pulse. The characteristic modes undergo decoupled 
single-mode squeezing transformations, and therefore they form a useful basis to describe the 
evolution of the entire multimode system. This provides an elegant and intuitive picture of quantum 
statistical properties of parametric fluorescence. We analyse the efficiency of detecting quadrature 
squeezing, and present results of numerical calculations for a realistic nonlinear medium. 
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I. INTRODUCTION 



Quadrature squeezing [l| is one of the most accessible optical processes that take us be- 
yond the semiclassical class of states of electromagnetic radiation, comprising coherent states 
and their statistical mixtures. The most elementary theoretical description of squeezing in- 
volves just a single radiation mode undergoing evolution governed by the unitary operator: 



where a is the annihilation operator of the mode, and ( is the interaction strength. Such evo- 
lution occurs in the process of parametric down-conversion, which can be realized by means 
of x nonlinearity in optical media. In this case, the constant £ involves the amplitude of 
the pump beam, the strength of the nonlinear coupling, and the length of the medium. A 
common experimental configuration is an optical parametric amplifier, where a short pump 
pulse passes once through the medium. This induces squeezing in a mode that is either 
initially left empty, or seeded by a coherent beam at half the pump frequency. Obviously, 
Eq. Q is a very crude approximation to the actual evolution of the system. The realis- 
tic down-conversion is in fact a subtle interplay of the spatio-temporal shape of the pump 
pulse with the phase-matching characteristics of the medium for different sets of frequencies 
and wave vectors. Nevertheless, the basic features of down-converted light predicted by the 
single-mode model explain well experimental observations. 

The purpose of this paper is to go beyond the single-mode approximation and to discuss 
theoretically the squeezing process in a more realistic model, which takes into account the 
above complexities. This results in multimode calculations, involving continuous degrees of 
freedom such as frequency and wave vectors, which can be carried out only by numerical 
means. Nevertheless, elementary single-mode squeezing does emerge even in the multimode 
model. This is a general consequence of the linearity of the evolution in the field operators 
and the associated existence of the Bloch-Messiah reduction @, S Q • However, the modes 
whose evolution can be described in such a simple way are of a very special form. We shall 
call them characteristic squeezing modes and analyse their properties. The shape of these 
modes, as well as the strength of the produced squeezing depend in a nontrivial way on 
the parameters of the system, including the amplitude and the shape of the pump pulse as 
well as the characteristics of the nonlinear medium. The description of realistic squeezing 
in terms of characteristic modes provides a simple intuitive picture of the process as well as 
means to optimize the amount of detected squeezing. 

This paper is organized as follows. First, in Sec. II we give the theoretical construction 
of the characteristic squeezing modes, and motivate it by the well-known Schmidt decompo- 
sition of a two-photon wave function. In Sec. Ill we perform the decomposition for a model 
Gaussian evolution of an optical parametric amplifier, and discuss its consequences for de- 
tecting quadrature squeezing by means of homodyne detection. Sec. IV presents numerical 
results for a down-conversion process realized in a nonlinear waveguide, and points out a 
simple scaling law for the squeezing parameters of the characteristic modes. Finally. Sec. V 
concludes the paper. 
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II. DECOMPOSITION 



We will assume here a one-dimensional evolution of the optical field, corresponding to 
propagation in a single-mode nonlinear waveguide. Let us first develop an intuition using 
the perturbative regime, extensively discussed in previous works in the context of spectral 
properties of biphotons jij E|. If the field at the down-converted frequencies is initially in 
the vacuum state |vac), the output state in the limit of weak pumping can be written as: 

U\vac) « | vac) + ^ J duj J dJ ^(u, u')a\uj)a ] (J)\v&c) . (2) 

The two-photon wave function ty(u,u)') is given by a product of of the spectral pump am- 
plitude at uj + oj' and the phase matching function at the down-converted frequencies. For 
a degenerate process, the wave function is symmetric with respect to the permutation of its 
arguments. Consequently, its Schmidt decomposition |7[ can be written in the form: 

oo 

^l{uj,J)=Y,^n^Wn^') (3) 
n=0 

where ( n are non- negative Schmidt coefficients, while are mutually orthogonal Schmidt 

functions, normalized to one. We have used here complex-conjugated Schmidt functions 
tpn( u ) i n or der to achieve correspondence with the general case discussed below. The Schmidt 
decomposition performed above suggests an introduction of another set of field operators 
defined by: 

b n = duj ijj n (u)a(uj) (4) 



These operators satisfy standard bosonic commutation relations: 

[b m , b n ] = 0, [b m , &£] = S mn (5) 

which follows directly from the fact that the Schmidt functions ipn{^) are orthonormal. The 
final wave function can now be written as: 

- oo 

U\vac) « |vac) + - ^ ( n {bl) 2 \v&c). (6) 

n=0 

It is seen that this wave function could be obtained from the perturbative expansion of the 
multimode squeezing operator of the form: 

Cn r /£t\2 12 



U = <$exp 2 ~bl] . (7) 

As the modes b n are mutually orthogonal, the evolution is effectively decoupled in the modes 
defined by the Schmidt functions, and if at the end we detect only one of these modes, then 
the evolution can be effectively treated in the single-mode approach. 

It turns out that the above decomposition exists not only in the perturbative regime, but, 
in a generalized form, also for arbitrarily intense multimode evolution as long as it remains 
linear in the field operators 0, HI H| • The non-perturbative case is most conveniently handled 
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in the Heisenberg picture, in which the final field operators are expressed in terms of the 
input ones. If the linearity condition is satisfied, the output field operators a out (oo) can be in 
general written in terms of the input field operators a m (oo) as a Bogoliubov transformation: 



a out ( 



oo) = doo' (C{oo,oo')a in (oo') + S(oo,oo')[a in {oo')^) (8) 



where the integral kernels C(oo,u') and S '(lo,oo') are the Green functions of the equations 
of motion for a concrete arrangement. In the perturbative regime, C (00,00') = 5(oo — oo') 
and the function S (00,00') is equal (up to a certain phase factor) to the biphoton spectrum 
amplitude ty(oo,uo'). Of course, the output operators must satisfy the same relations as the 
input field operators. This imposes the following conditions on the kernels: 



du"[C(u,u")S(u',u") - C (to', to") S (00,00")} = 0, 

doo"[C(oo, oo")C*(to', oo") - S(oo, oo")S*(oo', oo")] = 5(oo - oo'). (9) 

Because of these constraints, the singular-value expansions of the Green functions are not 
independent, but can be written using a joint set of functions and parameters as: 



C(u,J) = y^COShCn0r, 
n=0 
oo 

S(oo,oo') = 5>nhCnCW>')- (10) 

n=0 

Here <f) n (oo) and ip n (oo) are two orthonormal sets of functions, and ( n are real parameters. 
This decomposition is known as the Bloch-Messiah reduction Analogously to Eq. 
we can now define discrete sets of field operators according to: 

doo (p n (oo)a m (u) 

bT = J doj^ n {oj)a^\uj). (11) 

The Bogoliubov transformation introduced in Eq. (jHJ) induces now the following transfor- 
mation of the discrete field operators: 



tout 



S£cOShCn+(WsinKn. ( 12 ) 



It is seen that the evolution is completely decoupled between modes characterized by different 
indices n. Eq. (fT2j) is analogous to that of single-mode squeezing, and it generalizes the 
perturbative result derived in Eq. (JBJ) to the regime of multiple pair generation. Let us 
note that in general the shape of the input modes can be different from the shape of the 
corresponding output mode. As shown in Ref. j8j, with a non-chirped pump pulse, the 
input and output squeezing eigenmodes obey a relation ip n (oo) = </>*(cu), i.e. are the exact 
time reverse of each other. 
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III. HOMODYNE DETECTION 



Let us first discuss consequences of the multimode character of squeezing for homodyne 
detection with the local oscillator field described by a certain mode function ip LO {u). In 
order to expose the essence of the argument, we will assume the integration kernels in the 
following Gaussian form: 



C(u,u') = S(u-u'), 
S(u,u') = 



2N 



7T(5A 



exp 



(u + u' — u p ) 2 {u — u') 2 



25 2 



2A 2 



(13) 



These expressions can be justified on the grounds of the first-order perturbation theory 
and for simplicity we have assumed that the phase factors resulting from linear dispersive 
propagation have been compensated. The formula for S(u,u') involves u p as the central 
frequency of the pump pulse, 6 and A describing the frequency correlations within the photon 
pair as well as the biphoton bandwidth, and iV 1 equal to the average total number 
of generated photons. In a typical arrangement i5 C A and 5, given by the bandwidth 
of the pump pulse, defines the frequency anticorrelations within a biphoton, while A is 
determined by the crystal dispersion across the down-converted spectrum and specifies the 
overall bandwidth of the parametric fluorescence. 

The Gaussian form allows us to find analytically [j| the singular values as 



C n « sinh ( n 



N 



coshr 



tanh n r 



(14) 



where r = In y/K/S, while the corresponding characteristic functions are given by: 

4> n (u) = 1p n (u) = 



r, 



2 n n\,/n 



H n (r s (u - u p /2)) exp 



T, 



U, 



s 'u-^ 



(15) 



where H n (x) denotes the nth Hermite polynomial. Thus in our model the characteristic 
modes have the familiar shape of Hermite functions describing the eigenfunctions of a har- 
monic oscillator, with their width given by the parameter r s \ 



2 

5A' 



(16) 



Let us now assume that the local oscillator has a Gaussian spectral profile characterized by 
a bandwidth 5io- 



lpLo(u) 



exp 



U p \ 2 



(17) 



The local oscillator mode can be decomposed in the orthonormal basis of the characteristic 
output modes as: 



4>Lo(u) = ^Mnlpnju). 



(18) 



n=0 
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With the choice of the spectral phase as in Eq. (JT7jl the local oscillator pulse is matched 
temporally with the squeezed light, and only even characteristic squeezing modes appear in 
the decomposition owing to the defined parity of the modes. Furthermore, the expansion 
coefficients in Eq. ()18|) can be evaluated analytically Iol |: 



f , , / x ,* / x v / (2m)!tanh m r' , s 

M„^«*H = ^f^= (19) 

where r' = ln^o/v^A), and all M 2m +i = 0. As the modes ip n (u) are mutually orthogonal, 
the quadrature fluctuations detected in the local oscillator mode V'lo(^) are given by a sum 
of independent contributions from all the occupied characteristic modes. In each of these 
modes, the mean square quadrature noise is phase dependent, with the maximum (Qn+) 
and the minimum (Q 2 n _) given by 

(QkH jexp(±2Cn), (20) 

where the vacuum state noise is normalized to 1/4. Because a pure squeezed state is a mini- 
mum uncertainty state, the product of the minimum and maximum quadrature uncertainties 
is equal to that of the vacuum. In the model Gaussian case only even characteristic modes 
contribute, and furthermore the coefficients M 2m are real and thus have identical phases. 
Consequently we can set the phase of the local oscillator to a value which combines contri- 
butions from exclusively either squeezed or antisqueezed quadratures. Thus the fluctuations 
in the squeezed {Q 2 _) and the antisqueezed (Q\) quadratures are given by: 

1 oo 

(Ql) = - 4 J2 M Le ±2(2 -- (21) 



We see that the squeezing is still present, but is not necessarily of a minimum-uncertainty 
character. Uncertainty minimization occurs only if all the ( n , for which M n do not vanish, 
are equal. The latter requirement is fulfilled when the local oscillator mode is identical to 
one of the characteristic modes. However, even if this is not the case, minimum uncertainty 
is approximated remarkably well as long as the dominant terms in the series {M2 m } pick up 
characteristic modes with comparable squeezing parameters (,2m- 

In order to quantify the last statement, we make use of the fact that a deviation of a 
squeezed ensemble from the minimum-uncertainty character can be ascribed to a non-unit 
quantum efficiency 

-lG(Ql)m + MQl) + MQl)-l 
V A(Ql)+4(Ql)-2 • 

of the communication channel along which it has been delivered. Interpreting Eq. ()21|) in 
this fashion, we find the efficiency of the squeezed ensemble in the local oscillator mode. 
This quantity is plotted in Fig. Q where it is seen that high-efficiency squeezing is obtained 
in a broad span of modes. The spectral width of squeezed modes is limited from below 
by the master laser linewidth and from above by the phase matching range of the down- 
conversion arrangement. By choosing a sufficiently short crystal, one can control the upper 
limit, allowing generation of squeezed vacuum in arbitrarily short pulsed modes. Our results 
permit an intuitive interpretation if down- conversion is considered in the time representation. 
The two photons in a pair can be born at any moment within the duration of the pump 
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FIG. 1: Quantum efficiency of detected squeezing as a function of the local oscillator linewidth 
5lo for r = In y A/5 = 3. High quantum efficiency is observed for all 5los between 5 and A, with 
7? > 99 % for -2 < r' < 2. 



pulse, but necessarily at the same time. Therefore, all the infinitely short pulsed modes 
within the pump pulse must contain an even number of photons and are squeezed, and so 
are all the linear combinations of these modes — that is, all the modes that are shorter than 
the pump pulse. The existence of the upper limit A in the squeezing spectrum is explained 
by the nonzero length of the down-conversion crystal. Due to nonlinear dispersion (which 
limits the down-conversion spectrum), single-photon wavepackets propagating through the 
crystal diverge in time by A -1 . If the local oscillator pulse is chosen too short, it may happen 
that one photon in a pair is registered within the local oscillator mode, but the other one 
arrives either too soon or too late, leading to a reduction in detection efficiency. 

Most existing experiments on pulsed squeezing obtain the local oscillator from the master 
laser whose second harmonic serves as a pump for the down-conversion. Assuming that the 
second harmonic bandwidth is limited by the spectrum of the pump pulse, this situation 
corresponds to 5lo = 8 and consequently r' = — r. From Fig. Q] we find that such a choice 
of the local oscillator linewidth lies on the border of the high-efficiency region and is not 
optimal. Fig. |2] displays the quantum efficiency associated with such a measurement as a 
function of r. For small r's, this quantity approaches unity, but for highly eccentric biphoton 
spectra it appears to tend to a constant value of 0.86. 



IV. NUMERICAL RESULTS 

In this section we will discuss the non-perturbative regime for a realistic example of the 
parametric process in a nonlinear waveguide. For this purpose we will solve numerically one- 
dimensional propagation, described by the set of differential equations for the annihilation 
operators of monochromatic modes a(u; z) parameterized with the propagation distance z 
in the medium: 

9a ^ ] Z) = ik(Lu)a(cu; z) + — l — / du/ ^'^ z E p {u' + w)a\u'\ z). (23) 
az L NL E J 

The first term on the right-hand side represents the linear propagation of the field in the 
medium, while the second one is responsible for the nonlinear interaction. In the above 
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FIG. 2: Quantum efficiency of a pulsed squeezed state for the local oscillator pulse bandwidth 
equal to that of the master laser, depicted as a function of the parameter r characterizing frequency 
correlations. 

expression, E p {ui) is the spectral amplitude of the pump field, and k(ou) and k p {oj) are 
respectively the wave vectors of the signal and the pump fields. We have rescaled the 
spectral amplitude of the pump field by E = f duE p (uj), which allows us to introduce a 
single parameter L^l, expressed in the units of length, that characterizes the strength of 
nonlinear interaction. In particular, Lnl is inversely proportional to the amplitude of the 
pump field. In the case of a monochromatic pump and no phase mismatch between the 
pump and the signal fields, Lnl is a distance over which the root mean square quadrature 
noise scales by a factor of e. 

We have performed calculations for a down- conversion process taking place in a beta- 
barium borate (BBO) crystal and converting pump wavelengths centered at 400 nm to 
degenerate signal and idler wavelengths around 800 nm. The pump pulse was assumed to 
have a Gaussian shape of duration t p , corresponding to the spectral amplitude: 



with a flat spectral phase at half-way through the waveguide length. As the propagation 
equation (|23j) is linear in the field operators, it can be solved numerically by replacing a(u; z) 
and a'(cu;z) by c-numbers and correspondingly their complex conjugates. The complete 
Green functions are composed of solutions of the propagation equation for monochromatic 
inputs. Details of the numerical procedure have been described in Ref. 0. In Fig. El we 
depict the absolute values of the Green functions \C(u,u')\ and \S(u,ui')\ for several ratios 
L/L NL , with the length of the crystal L = 1 mm, and the pump pulse duration r p = 26 fs. It 
is seen that the analytical Gaussian form of the Green functions discussed in the preceding 
section should be indeed capable of modelling qualitatively the numerical results. 

Let us begin by inspecting the spectral intensity profiles for the characteristic modes 
ip n (uj). The examples shown in Fig. 0] demonstrate that the dependence of the profiles 
on the nonlinear interaction length is rather weak. Furthermore, it is easy to recognize 
the similarity of the intensity profiles to those described by the Hermite functions given in 
Eq. (|15j) . Nevertheless small deviations from the Hermitian shapes, especially breaking the 
defined parity of the characteristic modes, can have a dramatic impact on the amount of 
squeezing detectable with a restricted class of local oscillator shapes. The reason for this 
is that when detecting the squeezed quadrature, such local oscillator pulses would pick up 




(24) 
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an admixture of antisqueezed quadratures from odd characteristic modes, which can mask 
the effect of squeezing. It is also difficult to set all the phases in the characteristic-mode 
decomposition to the same value which selects only squeezed quadratures. This issue has 
been discussed in detail in Ref. 

The squeezing experienced by the nth characteristic mode is given by the squeezing 
parameter ( n . We found numerically that up to L/L NL = 15 there exist a simple scaling 
law: namely, that the squeezing parameter of any squeezed mode is inversely proportional 
to nonlinear length 

Cn = (25) 

We have introduced here a proportionality constant A n which we shall call the squeezing 
length. In order to discuss its physical meaning, let us first recall that in the case of a 
monochromatic pump and no phase mismatch between the pump and the signal fields, we 
would expect the squeezing parameter to be equal to the ratio of the interaction length to 
the nonlinear length ( = L/Lnl- However in the general multimode case studied here there 
is an appreciable phase mismatch and the pump field is a polychromatic ultrashort pulse. 
Consequently the squeezing parameters cannot be simply estimated as L/ Lnl- Nevertheless, 
Eq. (ffijj) suggests that one can introduce a notion of the effective crystal length A n for every 
characteristic mode, a quantity that depends on the mode number as well as on the pump and 
the phase matching properties of the medium. We propose the following way of intuitively 
understanding A n : for a given input mode shape n (u;) its interaction with the pump pulse 
remains phase matched, thus interacting efficiently, only over a finite length of the crystal, 
approximately equal to A n . This length is shorter for modes of broad spectral content as 
their various frequency components dephase more quickly from each other and from the 
pump pulse owing to crystal dispersion. These effects are clearly seen in Fig. El where we 
plot A n as a function of the mode number n. Indeed, A n s are smaller for modes with higher 
n, whose bandwidth is broader. Let us point out here that the Gaussian approximation fails 
to predict accurately the squeezing parameters at any Lnl, as the plot of the analytical 
values given by Eq. (fTlj) would be a straight line in the semi-logarithmic scale of Fig. |SJ 
This is a consequence of the fact that even in the perturbative regime the biphoton spectral 
amplitude is not Gaussian, but typically involves a much more slowly decaying sine function. 
Finally, let us note that first several modes have very similar squeezing lengths A n and thus 
squeezing parameters ( n . As discussed in the preceding section, this gives us freedom in the 
shapes of the local oscillator pulse that are capable of detecting maximum squeezing. 



V. CONCLUSIONS 

The Bloch-Messiah reduction is an indispensable theoretical tool in understanding the in- 
tricacies of optical parametric amplification in the broadband regime, when all the monochro- 
matic modes are coupled to one another through nonlinear interaction with the ultrashort 
pump pulse. The characteristic modes are a simple and elegant way to describe completely 
the quantum statistical properties of parametric fluorescence. They also form a convenient 
modal decomposition to analyse homodyne detection of squeezing. 

Numerical simulations presented in this paper show that the perturbative approach begins 
to fail when the ratio L/ Lnl approaches one. This threshold value corresponds to quadrature 
squeezing exceeding the order of 1/e, which is equivalent to over 8 dB in noise reduction. 
This suggests that the perturbative approximation should be quite sufficient for describing 
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FIG. 3: Green functions |C(u;,u/)| and \S(co, u/)| for the example of a 1 mm long waveguide in 
a BBO crystal and the increasing ratio L/Ljvx = 0.1 (upper pair), L/Lnl = 1 (centre pair), 
and L/Lnl = 10 (lower pair). The pump pulse duration is equal to t p = 26 fs and its central 
wavelength is 400 nm. Contours are drawn at 0.2, 0.4, and 0.6 of the maximum value. 
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FIG. 4: Spectral intensity profiles of the characteristic modes exhibiting strongest squeezing for 
interaction the lengths L/Lnl = 0.1 (solid), L/Ljvx = 1 (dashed), and L/Ljvx = 10 (dotted). 
Note that the dashed and the solid lines virtually overlap in the graphs. 
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FIG. 5: The squeezing lengths A n = CuLnl as a function of the mode number n for the validity 
range of the inverse linear scaling law, < L/Lnl < 15. 

squeezing experiments performed so far. As long as we stay within the perturbative regime, 
the Gaussian model developed in Sec. II shows that there is a certain degree of flexibility 
in selecting the bandwidth of the local oscillator, with the optimal range limited by the 
parameters S and A characterizing frequency correlations within biphotons. However, when 
we utilize the whole available bandwidth of the master laser for second-harmonic generation, 
extracting the local oscillator from the same master beam puts us with 5lo close to 5, which 
is right on the edge of the optimal region. One possible remedy is to limit the bandwidth of 
the pump pulse by generating second harmonic in a longer crystal, in the dispersion-limited 
regime. This would lower the value of 5 and place the master laser badwidth Slo used for 
homodyning inside the favourable region. 
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